Freezing of dynamical exponents in low dimensional random media 
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A particle in a random potential with logarithmic correlations in dimensions d = 1,2 is shown 
to undergo a dynamical transition at Td yn > 0. In d — 1 exact results demonstrate that Tdyn = T c , 
the static glass transition temperature, and that the dynamical exponent changes from z(T) — 
2 + 2(T C /T) 2 at high temperature to z(T) = 4T C /T in the glass phase. The same formulae are 
argued to hold in d — 2. Dynamical freezing is also predicted in the 2D random gauge XY model 
and related systems. In d = 1 a mapping between dynamics and statics is unveiled and freezing 
involves barriers as well as valleys. Anomalous scaling occurs in the creep dynamics. 



The motion of topological defects by thermal activa- 
tion over pinning barriers determines the slow glassy dy- 
namics in numerous disordered systems, e.g. domain 
walls in dirty magnets, vortices in superconductors, dis- 
locations in pinned lattices jjj. While treating many in- 
teracting extended defects in the presence of disorder re- 
mains a major challenge, progress can be made on the 
simpler, already non trivial problem of a single point de- 
fect. A model of great interest is a particle diffusing 
in a random potential with log-correlations, i.e barriers 
growing logarithmically with scale. In 2D it precisely 
describes a single vortex in a XY spin model with gaus- 
sian random gauge disorder B, and is relevant to a host 
of related systems e.g. vacancies in pancake lattices of 
layered 3D superconductors dislocations in 2D lat- 
tices with smooth disorder electrons on helium ||. 
As a prototype model of diffusion in a complex phase 
space or in random media, it is of broader interest to 
relaxation in glasses |J , transport in solids Q popula- 
tion biology non hcrmitian quantum mechanics |J 
and vortex glass dynamical scaling fig] . Similar models 
where used to study the dynamical generation of broad 
(e.g. power law) distributions of relaxation times and its 
relation to aging fll|| . Its static limit also appears in In- 
teger Quantum Hall (QH) transition studi es fll2|| and its 
quantum extension in QH bilayer systems ]l3j]. 

Early studies of this model, in the context of tracer 
diffusion in 2d potential flows, used RG methods pertur- 
bative in the disorder strength a Jl4|-|l(|. The position 
x(t) of a particle (in any d) satisfies a Langevin equation: 



x{t) = -W(x(t))+r)(t) 



(1) 



(rj(t)ri(t')) = 2T5(t — t') being a thermal noise (X and 
(X) respectively denote disorder and thermal averages). 
It was found that if correlations grow logarithmically: 



A(x - x') = (V(x) - V{x')) 2 ~ 4crln \x 



(2) 



the diffusion is anomalous with a temperature T depen- 
dent "dynamical exponent" z given by llaJlq]: 



2 + 2 



dT 2 



(3) 



Strikingly, this result for z was conjectured |lq] to hold 
to all orders in g = a/dT 2 . This was confirmed in d = 1 
and d = 2 by general arguments and a three loop cal- 
culation ijlTj (which also found 0(g 3 ) corrections to z in 
d > 3). Further support came from excellent agreement 
with simulations |l^], performed in d = 2 for < g < 0.8, 
and exact results in d = 1 [[l0| for the velocity v ~ f z ^ 1 
at small applied force /, with z again as in ([}]). 

Recently, however, the statics of this model has been 
investigated in several works in d = 2 P Jl2|[T9| -pl[| and 
any d p2{ . As a result we now know that there is in fact 
a transition at T = T c = \fo~Jd to a strong disorder, low 
T glass phase. This glass phase appears non trivial, as 
dominated by a few states, reminiscent of replica symme- 
try breaking (RSB). It is related, approximately, to the 
Random Energy Model (REM) 1 23 1 and, more closely, to 
the directed polymer on a Cayley tree (DPCT) |24j]. It is 
thus an outstanding problem to investigate whether this 
equilibrium transition has a dynamical counterpart, and 
how it can be compatible with (|§|). 

In this Letter we solve the apparent paradox. In d = 1 
we demonstrate that there is indeed a dynamical transi- 
tion at the static T c and obtain the dynamical exponent 
z(T) at all T. This is achieved through exact results 
and a real space RG (RSRG) method. Simple arguments 
and bounds also show that a dynamical transition occurs 
in d = 2, with analogous behaviour of z(T). A similar 
dynamical freezing is also predicted in the 2d random 
gauge XY model, since the XY phase and its boundary 
are dominated by the single (i.e dilute) vortex limit [p|,p2[. 

We first give a simple argument indicating that the 
1/T 2 divergence in (||) cannot hold at low temperature. 
First, in any given sample (in finite d), the characteristic 
time t for a particle to escape a region of size L satisfies: 

hxt w B/T < [V max (L) - V mm (L)]/T = B max /T (4) 

i.e is given as T —> by the Arrhenius estimate 
where B is the energy barrier encountered by the par- 
ticle, obviously bounded by the difference B max be- 
tween the absolute maximum V max (L) and minimum 
V m i n (L) of the potential in the region. Second, we 
now know |E0 12 2l| , |22]| that these extremal values sat- 
isfy V m i n (L) = — 2yJ a / d(d In L — ^lnlni) + SV where 



1 



7 is a universal number and SV has O(l) sample to sam- 
ple fluctuations. Thus, defining z = hit/ InL at large L 
p5[ this yields, in any dimension d, the bound: 



< 4Vad/T 



0. 



(5) 



Since (||) is exact perturbatively to all orders in g = a/T 2 
in d = 1, 2, implies that a dynamical transition m«st 
occur in d = 1, 2 at a finite temperature Td yn {d) > 0. 

In d = 1, which we now study, it is natural to guess 
that the upper bound in (||) gives the exact z(T) for all 
T < T c . Indeed one expects that the Arrhenius law (Q) 
should hold in the energy dominated glass phase and fur- 
thermore B should equal its upper bound in (|4|) since in 
d = 1 there is only a single path. We now confirm this 
with analytical results at all T using the first passage 
time approach |2(| . We also show that (i) in d = 1 there 
is a direct correspondence between dynamical (e.g. expo- 
nents) and static quantities (ii) the dynamical exponent 
z can be unambiguously defined. Our conclusions being 
independent of boundary conditions we choose them re- 
flecting at x = |^6| and define t as the first passage time 
at site x = L of a particle which starts at x = 0. Some of 
our results are valid for any potential landscape but are 
mostly applied to gaussian log-correlated potentials @ . 

As shown below it is sufficient to focus, in any given 
disorder realization V(x), on the (thermal) mean first 
passage time ti({V}) = (t) = r. We first obtain its 
typical behaviour (i.e in any sample), its large deviations 
being studied later. The exact formula for t\ {{V}) is [E6| : 



h({V}) = ~J dyj^ dx8(y-x)e 



[V(v)-V(x)]/T 



(6) 



where 9(x) is the step function. Since the canonical 
partition function for the same disorder realization is 
Z L {{V}) = J^dx exp(-V(x)/T), (§) strongly resem- 
bles two copies of the statics, one with disorder +V(x) 
(dominated by minima at low T), the other with —V(x) 
(dominated by maxima of V at low T). Dynamics and 
statics are thus directly connected by: lnjT^dU}) + 
h({-V})}} = \nZ L {{V})+\nZ L ({-V}) and also 

hiZ L/2 ({-V})+]nZ L/2 ({V}) < IniThdV})} 

<lnZ L ({-V}) + lnZ L ({U}). (7) 

If the ratio of upper to lower bound converges to unity 
in the thermodynamic limit, they uniquely determine the 
long time behavior in terms of the statics. This is the case 
for gaussian potentials with correlations growing at the 
most as a power of the logarithm of the scale. If the corre- 
lations do not grow faster than logarithmically with scale, 
then lnZidy}) ~ InL. If, additionally, the intensive 
free energy / = F L ({V})/ InL = -T 'In Z L ({U}) / 'In L is 
self-averaging, the dynamical exponent z is obtained as: 



z({V})= lim In h ({V})/ InL 

L— »oo 



-2//T, 



(8) 



and is also self averaging. This is the case both for un- 
correlated and log-correlated gaussian potentials. For the 
former (||) gives the normal diffusion value z — 2. For 
the latter, using the known results for / |ljj2^] we obtain 
our main result for the dynamical exponent: 

z{T) = z A (T) = 2 (1 + a/T 2 ) (for T > T c ) 
z(T) = 4V^/T (for T < T c ). (9) 

Thus a dynamical transition, away from the "annealed" 
value za(T) given by (^|), occurs at the same tempera- 
ture T^yn = T c = i/a as the equilibrium transition. At 
T c freezing occurs in the thermal configurations which 
dominate (i), simultaneously around minima and max- 
ima of the potential (i.e in the two copies), thus in 
the effective barrier. Interestingly, this transition coin- 
cides with the onset of logarithmic corrections. Indeed 
one can also characterize typical finite size fluctuations. 
Using that F L = /(InL - ^j(T) ln(lnL)) + 5F where 
5F has O(l) sample to sample fluctuations we find 
that ti({V})/Tt yp has a well defined 0(1) limit distri- 
bution at large L, the typical mean first passage time 
being T typ = L z{ - T \\nL)- a ^ with a(T) = 2j(T)^/T, 
7(T) = for T > T c but 7 (T C ) = \ and j(T) = § for 
T < T c . For faster growing (e.g. power law) correla- 
tions, the ratio between the bounds in Eq. (]?]) does not 
converge to one, and even the leading order of In t\ ({V}) 
still fluctuates at large L |26| as in the unbiased Sinai 
model (27). 

We now show as promised that z can be defined unam- 
biguously from the mean first passage time alone. In prin- 
ciple one could define a full set of dynamical exponents 
z p ({V}) = limi^oo lnt p ({V})/p\nL from higher ther- 
mal moments t p ({V}) = (t p ). Using the expression p6[: 



t P ({V}) = jr p 



-l) 



0<Xi, Vi <L i=1 



and 9 x ^ y = 9(x — y) < 1 we obtain, comparing with (jg), 
\nt p ({V})/p\nL < lnti({U})/lnL + lnp!/plnL, which 
together with the general inequality (t) p < (t p ) leads to: 



* P ({V})=zi({V}) 



and 



(10) 



for any integer p > 1 . This is because the thermal distri- 
bution of t has exponential decay and all moments p > 1 
are controlled at large L by the largest relaxation time in 
a given sample. Thus z — z\ characterizes the dynamics. 

By contrast, the distribution of escape times Pl{t) 
with disorder realization has broad tails extending in the 
region r ^> T typ ~ L Z ^ T \ where we obtained the estimate: 



Pl{t)*t« ^(^TeM-^^f 1 ). (ID 

Here fi = T/T c and (|ll|) is valid for (i) T < T c and (ii) 
for T > T c and z = In r/ InL > 4. For T > T c and 



2 



z(T) < z < 4 one has simply P l (t) « T ty 1 p (T typ /T) 1+p2 . 
This corresponds to a quadratic (plus linear) multifractal 
spectrum for rare occurrences of z. (|ll|) can be obtained 
by a Kosterlitz RG analysis of @, as in p3]. The result 
is a non linear Kolmogorov equation for Pj, (r) as a func- 
tion of InL, identical to the one describing the partition 
sum of two directed polymers on a Cayley tree, seeing 
opposite disorder, with constrained endpoints x < y. It 
yields ( pd[ ) up to log corrections, neglected here. More 
empirically (Jll|) can be obtained from the moments: 



I r L p , * 

(t) P =tp II [<lyi dx i e ^ ~ x i)\ ex P { \- A( - yi " x ^ 
"' t=l i,j'=l 

- ~A( W - tfJ ) - ^Afc - x J )]/2T 2 }. (12) 

which reads as a partition sum for 2p particles, p of type 
y (representing hills in the potential landscape) and p of 
type x (representing valleys) . Same type particles attract 
via a potential A(r)/2T while those of opposite type re- 
pel via — A(r)/T. This is similar to estimating Zp in the 
statics, except that there only one kind of particles (rep- 
resenting valleys) appears. Here hills and valleys play 
symmetric roles. ( |l2|) can be estimated for log-correlated 
potentials and for integer p > 1 as follows. At high T an 
"entropic" variational saddle point dominates (with all 
y's and x's far away O(L) from each other). At low T 
an "energetic" saddle point dominates (all y's close to- 
gether within O(l), all x's close together, y's and x's far 
away) . This yields the large L behaviour (universal since 
unaffected by changes in A(r) for small r): 



(t) p ~ L 2 p( 1+ -/ t2 ) for T>T c , p , 



L 2p(l/p+pa/T 2 ) 



for T <T r 



(13) 



i.e, as for Z p in the statics (of this model |19| , |l2| , 
the REM and the DPCT) there is a sequence of tran- 
sition temperatures T CiP = ^JpT c for the moments. 
One can check via a saddle point calculation of (t) p = 
J dTT p P L (T) that (H) is consistent with (0) and that 
the T CjP correspond to a change of behaviour from rare 
events to typical events dominance as the sad dle point 
crosses z sp — 4. By analytically continuing (t) p as p — > 
one can recover (^]) and the transition at T c in typical 
behaviour. The "entropic" saddle point still dominates 
for T > T c , but the "energetic" saddle point is replaced 
at T < T c by a one-step RSB ansatz. Each kind of "par- 
ticle" is arranged in p/m groups of m particles close by in 
space, while different groups are far away, with m = T '/T c 
at the optimum, extending the stati c PUp^ and the 
REM and DPCT replica calculations [|23|]. Compared to 
conventional static RSB which only involves "valleys" the 
interesting feature here is a nonequilibrium RSB which 
also involves "hills": indeed, near degeneracies of dis- 
tant barriers result, in the glass phase, in a nonequilib- 
rium splitting of the thermal distribution of the diffusing 



particle into a few packets in a single environment [E8| . 
These features are absent for weaker correlations (high T 
"entropic" saddle) or stronger ones p6| , where an "ener- 
getic" saddle without RSB dominates (degeneracies are 
subdominant in the Sinai landscape except in the pres- 
ence of a bias p7j). 

If an external force is applied, the creep velocity v re- 
lates to the disorder averaged mean escape time of re- 
gions of size Lq ~ 1// and is thus controlled by the 
annealed exponent za(T), distinct from z(T) at low T, 
a striking breakdown of naive dynamical scaling. Freez- 
ing manifests itself in large finite size corrections to v 
due to undersampling of the disorder average. From 
v -1 = ~(t) holds for fL 1, where (t) is given by 
the tilted landscape V(x) — fx. It can be approximated 
as the average over N samples of sizes L of the escape 
time in each sample, each distributed with Pj, (r) . A 
saddle point estimate shows that below T c , v ~ f ZA ~ Y 
for y = — lni/ In / > y* = 2/ fi 2 — 1, but that for smaller 
sizes 1 < y < y*, the typical v ~ f z m-i wnere th e expo- 
nent z m = v(l + y)/2 — y + 1 smoothly interpolates 
between the annealed and quenched one. 

The RSRG method, previously devised to describe dif- 
fusion in the Sinai landscape (for details see [p7| ) and in 
a broader class (29), allows to obtain complementary in- 
formation, e.g. about distribution of positions. Here, in 
the log-correlated landscape, it is implemented numeri- 
cally. From the original set of (i) alternating local ex- 
trema V(xi), (ii) their energy differences ("barriers" of 
heights Fi — \V{xi) — V(x i+ i)\) and (iii) the segments be- 
tween them ("bonds" of lengths ii = Xi+\ — Xi), one con- 
structs iteratively the "renormalized landscape at scale 
T" by removing as T increases all barriers between F and 
F + dT and merging the corresponding bonds. This dec- 
imation retains only the large barriers and deep valleys. 
The Arrhenius dynamics of a particle starting at xq is 
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FIG. 1. Probability distributions for a = |: (i) without 
and (ii) with additional short range disorder. Upper right: 
Prob r ((F - T)/F typ ) for 2 23 sites (i) with 2700 samples, ini- 
tial (r = 0, triangles) and asymptotic (r = 12, squares) (ii) 
with 1000 samples initial (r = 0, diamonds) and asymptotic 
(r = 18, circles). Lower left: asymptotic q(X = x/lr)\ (i) 
(r = 12, squares) (ii) (r = 18, circles). 
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FIG. 2. On same graph: dependence of (i) maximum bar- 
rier B m ax = Vmax — Vmin on system size L (ii) minimum 
barrier B — V left at decimation scale P on average bond 
length £r- Vertical axis: 4 — B/(i/a\nR), horizontal axis: 
ln(47r In R)/ In R with (i) B = B max , R = L (circles) and (ii) 
B = T, 7? = I r , L = 2 23 (squares) and L = 2 21 (triangles). A 
graph going to the origin means B ~ 4^/a\nR for i? — *• oo. 
Finite size corrections _B = 4- v /o : (lni? — i7ln(47rln_R)) with 
7 = 3/2 (dashed line) 7 = 1 (full line) are shown. 



any d > 2. Indeed, in order to escape the particle now 
only needs to find a set of saddles which connects to the 
boundary. A percolation and counting argument shows 
that it can do so by remaining within sites such that 
V(x)/lnL < 0. Thus the relevant barrier is bounded 
as B < —Vmin. In d — 2 this bound is likely to be 
saturated since the particle still finds deepest minima, 
yielding z = 2\f2ujT = 4T C /T. Since this expression 
matches (|j|) at the static T c in d = 2, a likely scenario 
is that Tdyn = T c and that the expression holds for all 
T < T c j3(J. A single vortex in a random gauge XY 
model will experience a similar dynamical freezing. 

To conclude we demonstrated dynamical transitions 
in d = 1,2. In d = 1 we found anomalous scaling of 
the creep velocity, novel freezing phenomena involving 
barriers and a finite T generalization of Arrhcnius law 
r^e~ 2Fi / T . Extensions will appear elsewhere. 



approximated by putting it at time t at the bottom 
of the bond which contains xq in the renormalized land- 
scape at T = Tint. The errors are small if the distri- 
bution of renormalized barriers is broad compared to T 
(infinitely broad in Sinai). Since here barriers remain 
finite, the method should be exact only as T/T c — > 0, 
i.e when the thermal packet is concentrated in a single 
well. Perturbative corrections in T/T c can be computed 
by considering the small occupation probability of neigh- 
boring secondary wells. In practice, the method works 
surprisingly better and gives the exact z(T) up to T c . 

We have found that under renormalization the proba- 
bility distributions for rescaled variables reach fixed point 
forms, namely Prob r ((^-r)/i^ yp ) -> P* ((F -T)/ F^p) 
for rescaled barriers (Fig. fy) and Probr(^/fr) — > 
Q*(i/l r ). Here Fr yp flows to a constant F^p w 4.4^ 
and £r — L/Nr is the average bond length. Since two 
barriers are decimated at each step the number of re- 
maining barriers Nr satisfies drNr = —2a-pN-p where 
a r = Prob r (F = T) -> a* = P*(0)/F typ , a constant. 
Thus the bond length grows as It ~ exp(2a*T) ~ t x l z 
and using T = Tint one recovers the dynamical expo- 
nent z = l/(2a*T). Numerically we find l/(2a r ) ~ 
4 v / cF(l — 0.25 ln(47r ln^r)/ ln£r) and thus a value of z con- 
sistent with @. The diffusion front, computed as in [p7| , 

converges to a scaling form as Probr(xt\00) — > £ r q(X = 
x/Ir), represented in Fig. |l|. Ir is thus the only relevant 
lengthscale, all moments of the displacement scaling as 



{x(tf) 



t k/z 



Finally we obtained good numerical 
evidence (see Fig. p that all above asymptotic scaling 
functions, as well as ceo = \foot* and F = F typ /^/a, 
do not change upon adding short-range disorder and are 
thus universal {Ir does change by a constant factor.) 

We now address d > 2. First we note that the bound 
(§) can be improved to z{T) < 2y/da/T as T for 
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